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^ , Abstract. We study the viscosity driven bar mode (Jacobi-like) instability of rapidly rotating quark matter stars 

O ' (strange stars) in general relativity. A triaxial, "bar shaped" compact star could be an efficient source of continuous 

' wave gravitational radiation in the frequency range of the forthcoming interferometric detectors. We locate the 

secular instability point along several constant baryon mass sequences of uniformly rotating strange stars described 
by the MIT bag model. Contrary to neutron stars, strange stars with r/|VK| (the ratio of the rotational kinetic 
energy to the absolute value of the gravitational potential energy) much lower than the corresponding value for 
^ . the mass-shed limit can be secularly unstable to bar mode formation if shear viscosity is high enough to damp 

QQ ' out any deviation from uniform rotation. The instability develops for a broad range of gravitational masses and 

, rotational frequencies of strange quark stars. It imposes strong constraints on the lower limit of the frequency at 

the innermost stable circular orbit around rapidly rotating strange stars. The above results are robust for all linear 
self-bound equations of state assuming the growth time of the instability is faster than the damping timescale. 
Whether the instability can grow or not depends on many different physical quantities (e.g. value of viscosities 
CO ■ (rather uncertain)). We discuss astrophysical scenarios where triaxial instabilities (r-mode and viscosity driven 

instability) could be relevant (a new born star, an old star spinning up by accretion) in strange stars described by 
the standard MIT bag model of normal quark matter. The spin evolution of a strange star strongly depends on 
^ l' the strange quark mass. Taking into account actual values of viscosities in strange quark matter and neglecting 

the magnetic field we show that Jacobi-like instability cannot develop in any astrophysicaly relevant temperature 
windows. The main result is that strange quark stars described by the MIT bag model can be accelerated to very 
high frequency in Low Mass X-ray binaries if the strange quark mass is consistent with values based on particle 
data tables. 
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>y , 1. Introduction et al. 1996, 1998). The m = 2, I = 2 "bar mode" insta- 
bility is the fastest-growing instability if a rotation rate is 

The first generation laser interferometric gravitational sufRcientlv high 
wave detectors such as LIGO, VIRGO and GEO600 

should be fully operationally soon and second genera- Neutron stars can be forced to rotate rigidly if viscosity 
tion detectors like LIGO II within a few years. Triaxial high enough to damp out any deviation from uniform 
instabilities of rotating compact stars (neutron stars, rotation. The viscosity-driven bar mode (Jacobi-hke) in- 
strange quark stars) can play an important role as emis- stability of neutron stars has been studied in general rela- 
sion mechanisms of gravitational waves in the frequency ^^^^^V Bonazzola, Frieben & Gourgoulhon (1996, 1998), 
range of these detectors (see e.g. Stergioulas 1998, 2003 or ^^o find that only the stiff-est equations of state (EOS) al- 
Andersson 2003 for a review). A rapidly rotating relativis- symmetry breaking. Since for "realistic" EOS of neu- 
tic star can spontaneously break its axial symmetry if the ^ron stars the ratio o{T/\W\ has a low value (usually lower 
ratio of the rotational kinetic energy to the absolute value ^^^"^ 0-14), this instabihty takes place only when a star is 
of the gravitational potential energy T/\W\ exceeds some rotating at a frequency close to the mass-shed limit and 
critical value (Gondek-Rosihska & Gourgoulhon 2002, Di ^^^^ high mass. 

Girolamo & Vietri 2002, Shapiro & Zane 1998, Bonazzola Strange quark stars (SS) are currently considered as a 

possible alternative to neutron stars as compact objects 

Send offprint requests to: D. Gondek-Rosiriska, e-mail: (see e.g. Weber 1999 for a review) . Claims of discovery 

Dorota.Gondek@obspm.fr have been recently announced (Drake et al. 2002), but the 
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evidence is not yet decisive (Walter & Lattimer 2002). It 
has been suggested by Gourgoulhon et al. (1999), Gondek- 
Rosinska ct al. (2GG0a, b; 2001a) that a triaxial instability 
could develop more easily in a uniformly rotating strange 
star than in a neutron star. It was shown that indepen- 
dent of any details of the EOS of strange matter, the ra- 
tio for self-bound objects like strange stars can be 
even more than twice as high as for any models of neutron 
stars. No stability analysis in general relativity has been 
performed for strange stars untill now, because numeri- 
cal constraints prevented investigators from treating stars 
with strongly discontinuous density profiles at the sur- 
face. In the present paper we perform a numerical study 
of the secular triaxial instability of rigidly rotating strange 
stars in general relativity. We use the code developed re- 
cently by Gondek-Rosiiiska & Gourgoulhon (2002) for lo- 
cating the secular bar mode instability points of uniformly 
rotating incompressible fluid bodies in general relativity. 
Numerous tests have been performed to assess the valid- 
ity of the method and the accuracy of the numerical code 
(Gondek-Rosihska & Gourgoulhon 2002). 

In this paper we shall focus on the following questions: 
for which gravitational masses and rotational frequencies 
of strange stars does the Jacobi-likc instability occur? 
Does gravitational radiation set strong limits on spin rate 
of strange stars and on the frequency at the innermost sta- 
ble circular orbit around rapidly rotating strange stars? 
What astrophysical scenarios does the viscosity driven in- 
stability allow one to develop? How do our conclusions 
depend on physical parameters of the MIT bag model? 

The plan of the paper is as follows. The description of 
previous studies of the bar mode instabilities is presented 
in Sec. 2. Section 3 and 4 are devoted to the description 
of the strange matter EOS used by us and of the numer- 
ical technique we employed, respectively. The numerical 
results and discussion are presented in two sections Sec. 
5 and 6. In Sec. 5 we examine the necessary condition 
for the onset of the instability: that the frequency of a 
corotating mode goes through zero in a frame co-rotating 
with the star. This condition yields the critical angular 
velocity (at a given compactness) for which the instability 
could start to operate, provided it is not damped by other 
mechanisms. The results arc robust for any linear self- 
bound EOS. In Sec 6. we discuss the astrophysical scenar- 
ios when triaxial instabilities could be relevant. This part 
includes an estimate of different timescales (uncertain). 
Finally, Sec. 7 provides a summary of our work. 

2. Previous studies of the bar mode instability 

Rapidly rotating compact stars may be subject to different 
kinds of bar mode non-axisymmetric rotational instabili- 
ties. 

There exist two different classes and corresponding 
time-scales for non-axisymmetric instabilities: a dynam- 
ical instability, growing on the hydrodynamical time- 
scale and secular instabilities which grow on much 
longer dissipation time-scales. In Newtonian theory, a 



rapidly rotating incompressible fluid body (a Maclaurin 
spheroid) is secularly unstable to a bar-mode formation, if 
T/|VF| > 0.1375 (Jacobi/Dedekind bifurcation point) and 
dynamically unstable if T/|iy| > 0.2738 (Chandrasekhar 
1969). A secular instability can grow only in the pres- 
ence of dissipativc mechanisms, such as shear viscos- 
ity (Roberts & Stewartson 1963) or gravitational ra- 
diation (the Chandrasekhar-Priedman-Schutz instability 
(hereafter CFS instability; Chandrasekhar 1970; Friedman 
& Schutz 1978; Friedman 1978). Viscosity dissipates me- 
chanical energy but preserves angular momentum. As 
shown by Christodoulou et al. (1995), the viscosity driven 
(Jacobi-like) bar mode instability appears only if the fluid 
circulation is not conserved. If on the contrary, the circu- 
lation is conserved (as in inviscid fluids submitted only to 
potential forces), but not the angular momentum, it is the 
gravitational radiation (Dedekind-like) instability which 
develops instead. If both viscosity and gravitational radi- 
ation are operative, they act against each other. Which 
of the two dissipative mechanisms dominates is highly de- 
pendent on the viscosities and the temperature of a star. 

The few numerical studies of secular instabilities in 
rapidly rotating stars in general relativity show that non- 
axisymmetric modes driven unstable by viscosity no longer 
coincide with those made unstable by gravitational radi- 
ation. As was shown by recent post-Newtonian (Shapiro 
& Zane 1998, Di Girolamo & Vietri 2002) and relativis- 
tic (Gondek-Rosihska & Gourgoulhon 2002) studies of 
viscosity-driven instability in rigidly rotating incompress- 
ible fluid bodies, general relativity weakens this bar mode 
instability: the values of r/|VF|, eccentricity and rota- 
tion rate increase at the onset of instability as the com- 
paction parameter GM/(i?circC^) increases. This stabiliz- 
ing tendency of general relativity is in agreement with 
previous studies of Jacobi-like instability of rigidly rotat- 
ing compressible (neutron) stars (Bonazzola, Friebcn & 
Gourgoulhon 1996, 1998). The stabilizing effect of relativ- 
ity is not very strong, the critical value of r/|VF| for rela- 
tivistic neutron star models is 0.14 (the Newtonian value 
is 0.1275 (Bonazzola et al. 1996) for compressible bodies) 
and 0.18 for very relativistic (with proper compaction 
parameter M/R = 0.25) incompressible fluid body. 

On the other hand, the gravitational-radiation driven 
instability is strengthened by relativity. As was first 
proved by Friedman & Schutz (1978) the Dedekind-like in- 
stability is generic in rotating stars - one can always find a 
mode that is unstable for a star with a given rotation rate. 
The most important pulsation modes for the CFS insta- 
bility are the r-mode and f-mode (Stergioulas 1998, 2003, 
Andcrsson 2003, Andersson & Kokkotas 2001, Stergioulas 
& Friedman 1998, Friedman and Lockitch 2002, for recent 
reviews of these instabilities) . For a wide range of realistic 
EOS, Morsink, Stergioulas & Blattnig (1999) find that the 
critical T/|iy| for the onset of the instability in the I = 2 
f-mode is only T/|VF| ~ 0.08 for models with mass 1.4M0 
and T/|iy| ^ 0.06 for maximum mass models. However, 
based on the relative strength of the f-mode and r-mode 
instabilities for NS, one can expect that the r-mode in- 
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stability will occur at lower rotation rates and in a larger 
temperature window than the f-mode instability, unless 
some damping mechanism favors a stronger damping at 
the r-mode frequency. The onset of the r-mode instability 
for SS has been studied by Madsen (2000) and Andersson, 
Jones & Kokkotas (2002). 

The dynamical instability against bar-mode deforma- 
tion of rapidly rotating polytropic stars was studied by 
Shibata, Baumgarte & Shapiro (2000) in full general rela- 
tivity and by Saijo et al. (2001) in the first post-Newtonian 
approximation. The authors considered differentially ro- 
tating proto- neutron stars. Their main result is that in 
general relativity the onset of dynamical instability oc- 
curs for a somewhat smaller critical value of (T/|W^|)crit ~ 
0.24 — 0.25 than the Newtonian value ~ 0.27 for incom- 
pressible Maclaurin spheroids. 

3. Equation of state 

The possibility of the existence of quark matter dates back 
to the early seventies. Bodmer (1971) remarked that mat- 
ter consisting of dcconfined up, down and strange quarks 
could be the absolute ground state of matter at zero pres- 
sure and temperature. If this is true then objects made of 
such matter, so-called "strange stars" , could exist (Wittcn 
1984). Typically, strange stars are modeled with an equa- 
tion of state based on the phenomenological MIT-bag 
model of quark matter (Alcock et al., 1986; Haensel et 
al., 1986), in which quark confinement is described by an 
energy term proportional to the volume (Farhi & Jaffe 
1984). 

In the framework of this model the quark matter is 
composed of massless u, d quarks, massive s quarks and 
electrons. There are three physical quantities entering the 
model: the mass of the strange quarks, Wg, the bag con- 
stant, B, and the strength of the QCD coupling constant 
a. 

The EOS of strange quark matter (SQM) is calculated 
using a perturbation expansion in a. The expansion is 
treated in the rcnormalization scheme. In particular, one 
has to use renormalized values of mg and a, which both 
depend on the selected value of the rcnormalization point, 
pj\ (sec, e.g., Farhi & Jaffe 1984). The value of pn is of 
the order of typical energies in the quark-system under 
consideration. Following Farhi & Jaffe (1984), we will use 
PR = 313 MeV. 

3.1. The strange quark mass 

The s quark mass given in the tables compiled by the 
Particle Group is the so called current quark mass, mi^-* = 
75—150 MeV (Groom et al. 2000). Such a mass can be used 
in the systems in which the energies of the relative quark- 
quark motion » 1 GeV, say 100 GeV, which is many 
orders of magnitude larger than typical energies in SQM. 
Therefore, the value of TOs has to renormalized using the 
rcnormalization group equation, and its value calculated 



at PR (Farhi & Jaffe 1984). The lowest-order (in a) rcnor- 
malization group equation for nis indicates that rcnormal- 
ization will significantly increase the s-quark mass, and 

(c) 

that with our choice of pR we will have rUs ~ 2ms • In 
what follows, we will use mgC^ = 200 MeV as a basic 
(standard) value of rris in SQM, but we will also present 
some results for TOs = 100 MeV, to visualize sensitivity of 
our results to the value of mg- 

3.2. Models of EOSs of strange matter 

Our basic EOS corresponds to the standard values of the 
MIT bag model parameters of SQM: 

- SQSl - the standard MIT bag model: TOgC^ = 
200 MeV, a = 0.2, B = 56 MeV/fm^; 

We shall also find the onset of the viscosity driven in- 
stability for two other MIT bag models to check how our 
results depend on model parameters: 

- SQSO - the simplified MIT bag model: nis = 0, a = 0; 
B = 60 MeV/fm^; 

- SQS2 - the "extreme" MIT bag model (relatively low 
strange quark mass and B but very high a) : mgC^ = 
100 MeV, a = 0.6, B = 40 MeV/fm^. 

4. Numerical codes 

4.1. Axisymmetric equilibrium stars 

We have calculated axisymmetric models of rotating 
strange stars using two different highly accurate relativis- 
tic codes. The first code developed by Gourgoulhon et al. 
(1999) is based on the multi-domain spectral method of 
Bonazzola et al. (1998). This code has been used previ- 
ously for calculating properties of rapidly rotating strange 
stars described by the MIT bag model (Gourgoulhon et 
al. 1999, Gondek-Rosihska et al. 2000a, 2001a, Zdunik et 
al. 2000a) and by the Dey et al. (1998) model of strange 
matter (Gondek-Rosihska et al. 2000b, 2001b). To com- 
pute very rapidly rotating axisymmetric stationary con- 
figurations (when the ratio of polar to equatorial radius is 
lower than 0.3 — 0.5 depending on the compaction param- 
eter - see Sec. III. A of Gondek-Rosihska & Gourgoulhon 
2002 for a discussion), where the first code has difficulties 
converging to a solution, we employ a second numerical 
code (Stergioulas & Friedman 1995, see Stergioulas 1998 
for a description). In this code the equilibrium models 
are obtained following the KEH method (Komatsu et al., 
1989), in which the field equations are converted to inte- 
gral equations using appropriate Green functions. Models 
computed by both codes agree very well in all computed 
properties (Stergioulas, Kluzniak & Bulik 1999). 

4.2. The secular bar mode instability 

Recently Gondek-Rosihska & Gourgoulhon (2002) used 
multi-domain method for calculating Jacobi-like instabil- 
ity of relativistic incompressible fluid bodies. They have 
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introduced surface-fitted coordinates, which enable them 
to treat the strongly discontinuous density profile at the 
surface of incompressible bodies. This avoids any Gibbs- 
like phenomenon and results in a very high precision, as 
demonstrated by comparison with the analytical result for 
the Newtonian Maclaurin-Jacobi bifurcation point, that 
their code has retrieved with a relative error of 10~^. Their 
results are in close agreement with the very recent post- 
Newtonian calculations of viscosity driven bar mode in- 
stability in rigidly rotating incompressible fluid bodies by 
Di Girolamo & Vietri (2002) (the relative difference is less 
than 

The method of finding the secular bar mode instabil- 
ity point presented there is an improvement on an ear- 
lier method by Bonazzola, Frieben & Gourgoulhon (1996, 
1998). The approach consists of perturbing an axisymmet- 
ric and stationary configuration and studying its evolution 
by constructing a scries of triaxial quasi-equilibrium con- 
figurations. If the perturbation is damped (resp. grows) 
during the interation, the equilibrium configuration is de- 
clared to be stable (resp. imstable) (sec Gondek-Rosinska 
& Gourgoulhon 2002 for a description). The above method 
does not require us to specify some value of viscosity. 
Whatever this value, the effect of viscosity is simulated 
by the rigid rotation profile that is imposed at each step 
in the iterative procedure. If one considers the iteration 
as mimicking some time evolution, this means that the 
time elapsed between two successive steps has been long 
enough for the actual viscosity to have rigidified the fluid 
flow. Consequently, using this method one cannot get the 
instability time scale, which depends upon the actual value 
of the viscosity, but only the instability. 

We use the same code as Gondek-Rosinska & 
Gourgoulhon (2002) for finding the viscosity driven in- 
stability points of rapidly rotating strange stars described 
by the MIT bag model. 

5. Calculations and results 

5.1. Effect of relativity on the secular bar mode 
instability 

In this section we are focusing on the effect of relativ- 
ity on the Jacobi-like instability in the case of strange 
quark stars. Strange stars in contrast to ordinary NS have 
a stiff EOS at the surface and are soft inside and their 
density profiles are flat. Rotating strange stars with small 
masses, lower than 0.1 Mq, can be described very well 
by a uniform-density fluid, i.e., the classical Maclaurin 
spheroids (see Amsterdamski et al. 2002 for comparison 
between fully relativistic calculations of rapidly rotating 
low mass strange stars and Maclaurin spheroids). Widely 
used indicators for the onset of viscosity driven instability 
are the ratio of the kinetic energy to the absolute value of 
the gravitational potential energy T/\W\ and the eccentric- 
ity e = (1 — (rp/req)^)^/^ where rp and req are the polar 
and equatorial coordinate radius. We have computed the 
T/|TV| according to the prescription of Friedman et al. 
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Fig. 1. Effect of general relativity on the critical value 
of the kinetic energy to the absolute value of the gravi- 
tational potential energy for different models of strange 
quark stars described by MIT bag model (solid line) and 
incompressible fluid stars calculated by Gondek-Rosinska 
& Gourgoulhon (2002) (dashed Hne). The ratio ofT/\W\ 
for compaction parameter GM/ (i?circC^) = corresponds 
to the classical Newtonian Jacobi/Dedekind bifurcation 
point, 0.1375. 
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Fig. 2. Critical values of stellar oblateness as a function 

of the proper compaction parameter for different models 
of strange quark stars described by MIT bag model (solid 
line) and incompressible fluid stars (dashed line) 



(1986): 

T/\W\ 



nj/2 



D.J/2 + Mp-M 



(1) 



where ft is the angular velocity, J the stellar angular 
momentum, M is the gravitational mass and Mp is the 
"proper " internal energy. A rotating configuration be- 
comes unstable to a bar-mode formation if the ratio of 
T/|W| or the eccentricity exceed some critical value. For 
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incompressible Newtonian bodies, the critical values of 
T/\W\ and e are 0.1375 and 0.8126 respectively. The 
Newtonian value of r/|M^| is 0.1275 (Bonazzola et al. 
1996) for compressible bodies. 

General relativity weakens the viscosity driven bar 
mode instability, but the stabilizing effect is not 
very strong (Gondek-Rosihska & Gourgoulhon 2002, Di 
Girolamo & Vietri 2002, Bonazzola et al. 1996, 1998). 
The ratio r/lVFl for relativistic neutron star models is 
~ 0.14 and for highly relativistic (GM/iRdrcC^) = 0.25) 
uniformly rotating incompressible bodies is only 30% 
higher than the Newtonian value (Gondek-Rosihska & 
Gourgoulhon 2002). It was found by Gondek-Rosihska & 
Gourgoulhon (2002) that for rigidly rotating homogeneous 
fluid bodies the dependence of (r/|W^|)crit on the proper 
compactness parameter x :— GM / (RdrcC^) can be very 
well approximated by the function 



iT/\W\% 



(T/|M^|)„.it,Ncwt + 0.148a;(a;-|-l) 



We locate the onset of instability (see section 4.2 
for description of our method) along several constant 
baryon mass sequences of uniformly rotating axisymmet- 
ric strange quark stars described by different MIT bag 
models (see section 3.) for the compaction parameter 
GM/(i?circC^) = - 0.25. To show the role of relativis- 
tic effects we plot in Fig. ^the ratio of r/|W| as a func- 
tion of the proper compaction parameter at the secular 
instability points for strange stars. All 3 models described 
in section 3 are represented by a solid line. We compare 
our results with those obtained by Gondek-Rosinska & 
Gourgoulhon (2002) for relativistic incompressible bodies 
(dashed hnc). We see the same dependence of the crit- 
ical value of (r/|T4^|)crit on compactness in both cases. 
In the Newtonian limit {CM/ {RdicC^) — 0) the ratio of 
{T/\W\Ut is 0.1375. 

In FigElwe plot the oblateness of the star (the ratio of 
fp/feq) versus the proper compaction parameter for both 
MIT bag model strange stars (solid line) and incompress- 
ible fluid bodies (dashed line). It was found by Gondek- 
Rosihska & Gourgoulhon (2002) that the critical value of 
the eccentricity depends very weakly on the degree of rela- 
tivity and for GM / [RcucC^) — 0.25 is only 2% larger than 
the Newtonian value at the onset for the secular bar mode 
instability. In the case of strange stars the dependence is 
a little bit stronger. For massive strange stars the density 
profile is no longer completely flat - for example the ratio 
of central density to surface density for a strange star is 
1.25 and 1.8 for proper compaction parameters 0.1 and 
0.2 respectively. 

5.2. Results for the standard MIT bag model 

Here we present results for the standard MIT bag model, 
SQSl, where the mass of strange quark is taken to be 
200 MeV/c^, the coupling constant 0.2 and the bag con- 
stant B = 56 MeV/fm"^. We show that the results pre- 
sented in this section can be applied to the other MIT 
bag models of SS to a very good approximation (Sec 5.3). 



5.2.1. Constant baryon mass sequences 

We construct constant baryon mass sequences of rotat- 
ing strange stars described by the SQSl model, i.e. the 
so-called evolutionary sequences. This is relevant for the 
evolution of an isolated strange star slowly losing its en- 
ergy and angular momentum via gravitational and electro- 
magnetic radiation. In general we cannot neglect accretion 
on a compact star (e.g. an accretion from supernova on a 
just born strange star or from a companion in low mass 
X-ray binaries) . 

We compute both normal (with Afb = 0.8, 1.6, 2.IM0) 
and supramassive sequences (2.5, 3.0 Mq). A sequence 
is called normal if it terminates at the zero angular mo- 
mentum limit with a static, spherically symmetric solu- 
tion, and it is called a supramassive sequence if it does 
not. The boundary between these two sequences is the se- 
quence with the maximum baryon mass of a static configu- 
ration, which is 2.I7M0 for the SQSl model. We consider 
configurations stable with respect to axisymmetric per- 
turbations, i.e for a supramassive constant baryon mass 
sequence a model is stable to quasi-radial perturbations, 

- — I < (or ( — — I < 0) where Mb is the 

baryon mass of the star and pc is the central density of 
the star (Friedman, Ipser & Sorkin 1988). 



if 



The dependence of angular momentum versus central 
density as well versus rotational frequency for the constant 
baryon sequences is shown in Fig.|31 Each sequence is rep- 
resented by a solid line and is labeled by the baryon mass 
as well as by the gravitational mass of the static configu- 
ration (in parentheses), if it exists. The gravitational mass 
of a strange star rotating at the Keplcrian limit is only a 
few percent (up to 7%) higher than M of the static mem- 
ber of the sequence (for the infiuence of rotation on the 
gravitational mass see Figs. 2 in Gondek-Rosihska et al. 
2000b and Gondek-Rosihska et al. 2001a). Marginally sta- 
ble configurations to quasi-radial perturbation are shown 
as a dot-dashed line. The angular momentum increases 
along each sequence from J = for static configurations, 
or Jmin on the dash-dotted line for supramassive stars, to 
"/max for the Keplcrian configurations (thin short dashed 
line, as well as the crosses for each considered sequence). 
The star in a sequence reaches the mass-shedding limit if 
the velocity at the equator of a rotating star is equal to 
the velocity of an orbiting particle at the surface of the 
star. It is worth noting (Fig. that for any given baryon 
mass, the maximal rotating configuration is not mass- 
shedding one. The mass-shed configurations are reached 
due to the increase of the equatorial radius relative to the 
deformation of the rotating star. Such a phenomenon was 
discussed by Zdunik et al. (2000b) in the case of normal 
sequences of MIT bag model strange stars and Gondek- 
Rosihska et al. 2000 for the Dey et. al (1998) EOS. This 
feature is characteristic of stars described by a self-bound 
linear EOS. If a strange star is born with high angular mo- 
mentum (close to the maximum available) it will spin up, 
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Fig. 3. Angular momentum as a function of central density (left) and the rotation frequency / = Q./{2'k) (right) 
along sequences of constant baryon mass (solid lines). Each solid line is labeled by this baryon mass in solar mass 
units, as well as (in parentheses) by the gravitational mass of the static configuration, if it exists. The dash-dotted 
line indicates stars marginally stable to axisymmetric perturbations. The Keplerian configurations are shown as thin 
short-dashed line (also as crosses for each of considered sequences). The long dashed line indicates stars marginally 
stable to non- axisymmetric perturbations. Thick solid lines represent the configurations secularly unstable to the bar 
mode formation taking into account viscosity as a leading dissipative mechanism 



losing the angular momentum (e.g by emitting of gravita- 
tional waves). 

Marginally stable configurations to non-axisymmetric 
perturbation are shown as a long-dashed line (as well as a 
filled circle for each considered constant baryon mass se- 
quence). As mentioned in Sec. 4, the multi-domain spectral 
technique with surface-fitted coordinates that we employ 
has some trouble when rp/req < 1/3 — 0.5 depending on 
the compaction parameter. To compute very oblate and 
highly relativistic (with high compaction parameter) ro- 
tating axisymmetric stationary configurations we use the 
Stergioulas & Friedman (1995) code, but we cannot lo- 
cate the Jacobi-like bifurcation point in this case. The 
thin part of the long dashed line (for / > 1.3 kHz for the 
SQSl model) represents marginally stable configurations 
found in the following way - we noticed that the depen- 
dence of central enthalpy on rotational frequency f{Hc 
and / are two input parameters needed for finding a rotat- 
ing star described by given equation of state) at the onset 
of instability is very well approximated by the following 
formula: = 0.47/[kHz] - 0.45 for SQSl model. We ex- 
trapolated this formula to higher rotational frequency in 
order to find central enthalpy to compute the axisymmet- 
ric configuration. 

Thick parts of the solid lines, between the filled circle 
and the cross, correspond to configurations secularly un- 
stable to the bar mode formation. These configurations 
can be good emitters of gravitational waves. However 
a strange quark star described by this model must ro- 
tate very rapidly, with frequency higher than 1.1 kHz for 
1.4 Mq, to allow symmetry breaking. If viscosity is high 
(see discussion in Sec 6), an accreting strange star (be- 



longing to a normal sequence) in low mass X-ray bina- 
ries (LMXBs) will never reach the mass-shedding limit. It 
can be accelerated to the critical frequency (which can be 
higher than the Keplerian frequency - see Fig. 3), corre- 
sponding to the marginally stable configuration to the bar 
mode formation instead. 

5.2.2. Masses and rotational frequencies for secularly 
unstable strange quark stars 

The crucial question from the observational point of view 
is - for which masses and frequencies the triaxial instability 
in the case of strange stars can develop. In order to an- 
swer this question we plot the value of T/|Vl^| versus grav- 
itational mass for several constant frequency / — ri/27r 
sequences (thin solid lines) of uniformly rotating strange 
stars in Figure 01 On one end of each sequence there is a 
mass-shedding configuration (with the lowest central den- 
sity in the sequence), and on the other end the last sta- 
ble configuration with respect to axisymmetric perturba- 
tions (the densest object in the sequence). The dotted line 
separates normal strange stars from supramassive strange 
stars. 

We see that the ratio of r/|iy| for strange quark stars 
rotating at the mass-shedding limit is very large and can 
be even higher than T/\W\ = 0.32. The T/\W\ is very 
high for all linear self-bound EOS (see Gondek-Rosinska 
et al., 2000b for different MIT bag models and Gondek- 
Rosinska et al., 2000a for strange stars described by the 
Dey model). The values of for strange stars are 

much larger than for any neutron star models (for which 
the maximum T/|iy| ranges typically from 0.1 to 0.14 - 
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Fig. 4. The ratio of the kinetic energy to the absolute 
value of the gravitational potential energy as a function 
of gravitational mass for the standard MIT bag model of 
strange quark stars. The Newtonian limit is obtained for 
M ^ {J ^ 0, Pc ^ Po) (a rotating strange star can 
be considered as the Maclaurin spheroid). The solid lines 
correspond to sequences of configurations rotating with 
constant rotational frequency. The rotational frequency is 
labeled close to each line. On one end of each sequence 
there is a mass-shed configuration (short-dashed line) and 
on the other end the last stable configuration with re- 
spect to axisymmetric perturbations (dash-dotted fine). 
The dotted line separates normal strange stars from supra- 
massive strange stars. The meaning of the other lines is 
the same as in Fig. 3. Configurations with T/|VK| higher 
than represented by long-dashed line (e.g. thick parts of 
solid lines) are secularly unstable to the bar mode forma- 
tion taking into account viscosity as a leading dissipative 
mechanism. 



see e.g. Cook et al., 1994). In contrast to neutron stars, the 
T/|iy| for strange stars increases with decreasing stellar 
mass for all models. Very high values of T/|VF| are charac- 
teristic not only of supramassive stars but also of rapidly 
rotating strange stars with small and moderate masses. 
The large value of T/|VF| results from a flat density pro- 
file combined with strong equatorial flattening of rapidly 
rotating strange stars. The long-dashed line represents the 
onset of the secular bar mode instability (critical values 
of and e). Conflgurations with T/IM^I higher than 

represented by this line are unstable with respect to non- 
axisymmetric perturbation and can be good emitters of 
gravitational waves (thick solid lines). 

From Figs 01 and El we see that strange stars could 
break their symmetry far away from the Keplerian limit 
and for broad ranges of masses and rotational frequencies. 
The critical value of T/\W\ is in the range 0.1375 - 0.18 
for very low mass rotating strange stars (in the Newtonian 
limit very well approximated by Maclaurin spheroids) and 
supramassive stars respectively. Those values are much 




11 13 

R [km] 

Fig. 5. Gravitational mass vs. radius for stars described 
by the standard MIT bag model. All lines have the same 
meaning as in Fig. 01 Models located to the right of the 
long dashed line are secularly unstable to the viscosity 
driven bar mode formation. 



lower than the Newtonian value of the dynamical instabil- 
ity against bar-mode deformation, 0.27 (according to our 
knowledge no stability analysis has been performed yet for 
uniformly rotating stars in general relativity). Existence 
of a crust only marginally affects our results. It has lit- 
tle influence on location of the viscosity driven instability 
points. According to Zdunik, Haensel and Gourgoulhon 
(2001) the maximal values of r/|iy| at mass-shedding 
limit are only 20% lower for rapidly rotating strange stars 
with a maximal crust than the corresponding values for 
bare strange stars described by the same MIT bag model. 
From Fig. 3 we see that a high rotational frequency, 
higher than 0.95 kHz (> 1.1 kHz for 1.4 Mq), is needed 
for a star described by the standard MIT bag model to 
reach the instability region. If we consider a model with 
TOg = 200 MeV/c^, but with high coupling constants 
a = 0.6 and low bag constants B = 36.4 MeV/fm'^ then 
the critical frequency for which the viscosity driven bar 
mode instability could occur is always > 0.77 kHz and is 
0.88 kHz for a canonical strange star. 

5.3. Results for other MIT bag models 

In Fig Elwe show gravitational masses and rotational fre- 
quency for which triaxial instabilities in rotating strange 
quark stars can develop. Each line corresponds to the on- 
set of instability for particular MIT bag model described 
in Sec 3. Strange stars rotating with a frequency higher 
than represented by this line are unstable with respect to 
non- axisymmetric perturbations and can be good emitters 
of gravitational waves. We see that viscosity driven insta- 
bility develops for a broad range of gravitational masses of 
strange stars (from that of a planetoid to about three so- 
lar masses) independently of the model. However, a high 
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Fig. 6. Rotational frequency vs. gravitational mass for 
secular ( viscosity driven) instability configurations de- 
scribed by different MIT bag models. Stars rotating with 
higher frequency than critical one (shown by each curve) 
are secularly unstable to non-axisymmetric perturbations. 
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Fig. 7. Rescaled Figure El to see how results depend on 
the a parameter at equation of state given by eq. (1). The 
onset points of stability for incompressible fluid bodies is 
shown for comparison (thin solid line). All other lines and 
symbols have the same meaning as in Fig. |S1 

rotational frequency (always higher than the frequency 
~ 0.64 kHz of the fastest observed pulsar) is needed for a 
star to reach the instability region. The higher the gravi- 
tational mass the higher is the frequency needed to form a 
"bar-shaped star". The most comfortable situation is for 
MIT bag models with low MIT bag constant (e.g. SQS2) 
and for small stellar masses. However one should keep 
in mind that detailed simulations of supernovae do not 
predict remnant masses less than 1.2Mq (Timmes et al. 
1996). 

It was shown that strange stars described by different 
models of strange matter can be very well approximated 



(with accuracy 1 — 2%) by a linear function (Zdunik 2000, 
Gondek-Rosihska et al. 2000a) 

P = a(p-po)c^ (1) 

where P is the pressure, p the mass-energy density, and 
c the speed of light. Both po E^nd a are functions of the 
physical constants B, nis and a for the MIT bag model. 
In general, eq. (1) corresponds to self-bound matter with 
density po at zero pressure and with a fixed sound veloc- 
ity (Vac) at all pressures. For the MIT bag model the 
density of strange matter at zero pressure is in the range 
~ 3 X 10" - 6.4 X 10" g cm-3 and a between 0.289 (for 
a strange quark mass of msC^~ 250 MeV) and 1/3 (for 
massless quarks). For a fixed value of a, all stellar pa- 
rameters are subject to scaling relations with appropriate 

1/2 

powers of po, / Pq , where / denotes either of the fre- 
quencies (rotational or orbital) , and M, i? cx pg (see 
e.g Witten 1984; Zdunik 2000). The scaling law for stel- 
lar parameters with a depends on rotational frequency 
(Gondek-Rosihska et al 2003), for approximate scalings 
with a at the limiting case of rotation at mass-shedding 
frequencies see Stergioulas, Kluzniak & Bulik (1999). 

In Fig. [3 the dependence of rotational frequency 
vs. mass for marginally stable strange stars to non- 
axisymmetric perturbations (Fig. has been scaled out 
with appropriate powers of po = 10^'' g cm~^. We see 
that the function f{M) for onset points of instability very 
weakly depends on a (a changes not very much in MIT 
bag model). The results obtained for one model can be 
used for other MIT bag models with the use of the scaling 
with po- The critical points for incompressible fluid stars 
with po — const = 10^^ [g cm~'^] are shown for comparison 
(thin solid line). 

From FigUlwe see that independently of any details of 
SS EOS, the lower the stellar mass the lower rotational fre- 
quency /crit needed to reach the instability region. Taking 
this into account we can find the absolute lower limit on 
rotational frequency of SS when viscosity driven instabil- 
ity can develop. In the Newtonian limit, when M — > 0, 
J ^ 0, pc ^ Po, a rotating strange star can be consid- 
ered as the Maclaurin spheroid for which the rotational 
frequency depends only on density and eccentricity (or 
r/|W^|) (Chandrasekhar 1969). The square of the angular 
velocity at the onset point of viscosity driven instability 
(e = 0.8126) is 



in other words for a strange star with gravitational mass 
M the frequency -Jm^ 0.4458 [kHz]. The lowest 
value of /crit is obtained for the lowest value of po (ful- 
filling the condition of stability neutrons to spontaneous 
fusion into droplets of u and d quarks as well as that en- 
ergy per baryon for strange quark matter must be below 
930.4 MeV). For the MIT bag model we can find that the 
lowest rotational frequency at which triaxial instability 
can develop is ~ 0.75 kHz. 
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Fig. 8. The secularly unstable configurations to the bar 
mode formation (thick solid lines between a heavy black 
circle and a cross) on the maximum orbital frequency ver- 
sus frequency of the rotation plane. The long dashed line 
indicates stars marginally stable to nonaxisymmetric per- 
turbations taking into account viscosity as a leading dis- 
sipative mechanism. One sequence for a very low mass 
star (with Mb = O.OIM©) is shown, the critical point 
on this sequence for Newtonian dynamical instability to 
non-axisymmetric perturbations is indicated by an aster- 
isk and dynamically unstable configurations are denoted 
with the dotted line. 

5.3.1. Lower limits on the maximum orbital frequency 
around rotating strange quark stars and 
kilohertz QPOs 

In this section we show that triaxial instability is im- 
portant for interpretation of maximal orbital frequency 
of rapidly rotating strange stars in LMXBs. We focus 
on the MIT bag model but our results are robust for 
other linear self-bound EOS. The spin-up of a relativistic 
star crucially depends on the maximum orbital frequency 
around it, as do a host of other high energy accretion 
phenomena in low mass X-ray binaries, including quasi- 
periodic oscillations (QPOs) in the X-ray flux. Kluzniak 
et al.(1990) suggested that the mass of a neutron star 
may be derived by observing the maximum orbital fre- 
quency, if it occurs in the marginally stable orbit, and 
that the low frequency QPOs occurring in X-ray pulsars 
will have their counterpart in LMXBs, at frequencies in 
the kHz range. Such kHz QPOs have indeed been discov- 
ered and their frequency used to derive mass values of 
about 2M0, under the stated assumption (Kaaret et al. 
1997, Zhang et al. 1998, Kluzniak 1998, Bulik 2000). In 
fact, the QPO frequency may correspond to a larger orbit, 
and hence a smaller mass (Kluzniak & Abramowicz 2002). 
The question of whether quark stars may have maximum 
orbital frequencies as low as 1.07 kHz in 4U 1820-30 (re- 
ported by Zhang et al. 1998) has also been investigated 
(Bulik et al. 1999a,b, Stergioulas et al. 1999, Zdunik et 



al. 2000a, b, Gondek-Rosihska et al. 2001a, b). The main 
result of these considerations was that the lowest orbital 
frequency at the innermost marginally stable orbit (ISCO) 
is attained either for non-rotating massive configurations 
close to their maximum mass limit, or for configurations at 
the Keplerian mass-shedding limit. In order to obtain an 
orbital frequency in the ISCO as low as 1.07 kHz for slowly 
rotating strange stars one has to consider a specific set of 
(rather extreme) MIT bag model parameters (Zdunik et 
al. 2000a, b) allowing gravitational masses higher than 2.2 
Mq. In contrast, for stars rotating close to the Keplerian 
limit, a low orbital frequency at the ISCO can be obtained 
for a broad range of stellar masses (Stergioulas et al., 1999, 
see Gondek-Rosihska et al. 2001a for discussion). 

We present the dependence of the maximum orbital 
frequency (the frequency at the marginally stable orbit 
if it exists or the Keplerian frequency at the stellar sur- 
face, /orb, max = /orb(-R*)) on the frequency of rotation f 
for the set of constant baryon mass sequences of strange 
stars in Fig. |H1 Marginally stable orbits exist for a broad 
range of masses and rotational rates (see e.g. Fig 1 in 
Gondek-Rosihska & Kluzniak 2002). As it turns out, the 
very low mass limit, represented in Figure |S1 by the case 
of Mb = O.OIM©, is described well by the Maclaurin 
spheroid for which analytic formulae for orbital frequen- 
cies have been derived and found to agree with the rela- 
tivistic numerical models (Amsterdamski et al., 2002). For 
this sequence, an ISCO exists only for rotation rates larger 
than ^ 950 Hz. We indicate the onset of the (Newtonian) 
dynamical instability to non-axisymmetric perturbations 
by an asterisk. In the Newtonian limit, the gap between 
the marginally stable orbit and the stellar surface is pro- 
duced by the oblatcncss of the rapidly rotating low-mass 
quark star (Kluzniak, Bulik & Gondek-Rosihska, 2001; 
Amsterdamski et al., 2002; Zdunik & Gourgoulhon 2001). 
We see that very rapidly rotating strange stars are sec- 
ularly unstable to bar mode formation if the viscosity is 
high enough. The lowest frequency in the marginally sta- 
ble orbit is obtained for slowly rotating massive strange 
stars or a marginally stable configuration to nonaxisym- 
metric perturbation (long dashed line). If the mass of 
the strange quark is relatively small, mgC^ ~ 100 MeV, 
then the r-modes instability could play an important role 
as an emission mechanism of gravitational waves and a 
mechanism spinning down rotation rates of strange stars 
(Andersson et. al. 2002). 

6. Rotating strange quark stars as sources of 
gravitational waves 

The interesting question is whether the viscosity-driven 
instability domain can be reached by any astrophysical 
process. In section 5 we located the onset of instability 
(when the frequency of an / = — m mode goes through 
zero in the corotating frame) along constant baryon se- 
quences of rotating strange quark stars. This necessary 
condition yields the critical angular velocity for a given 
gravitational mass of a star. From our numerical results, 
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the rotation frequency has to be above 0.8 — 1.2 kHz for 
strange quark stars described by the MIT bag model of 
normal quark matter (any kind of superconductivity being 
excluded). The instability can develop for a broad range 
of gravitational masses of rotating SS. 

The above results were obtained under the assumption 
of infinite shear viscosity. The necessary requirement is 
that the 'real' viscosity has to be high enough to dissipate 
energy and let the instability develop along the Jacobi 
path. 

As we mentioned in Sec 2 if both viscosity and grav- 
itational radiation are operative, they act against each 
other. Since the viscosity-driven mode is triaxial it tends 
to be damped by gravitational-wave emission, and since 
the gravitational-wave driven mode involves differential 
rotation, it is damped by viscosity. The instability can be 
important if the typical instability growth rate (t^ shear 
viscosity timescale) is faster than any damping mecha- 
nisms. 

It is important to emphasize that the spin evolution of 
strange stars depends on many different physical mecha- 
nism, concerning e.g. viscosities of dense matter, cooling 
evolution, the influence of magnetic field on instabilities 
etc, for which our understanding is still unsatisfactory. 
Below we consider different timescales (shear viscosity, 
bulk viscosity, cooling timescale, gravitational radiation) 
and provide a crude estimate of the stellar temperature 
for which viscous effects could dominate. Viscosity driven 
instability is divergence free so bulk viscosity does not in- 
fluence it. The shear and bulk viscosity (axial modes are 
divergence free but mode coupling should be taken into ac- 
count) of strange quark star matter is able to suppress the 
growth of the CFS-instability except when the star passes 
through a certain temperature window. The Z = m = 2 r- 
mode is expected to lead to the strongest CFS-instability 
so we will take into account only this mode in our discus- 
sion. 

6.1. Viscous timescales connected to sliear viscosity 

Shear viscosity of normal SQM has been calculated by 
Heiselberg and Pethick (1993): 

rj^5x 10^5 (^H^ (^A^ g cm-i , (2) 

where Tg = T/10^ K. Assuming p — 4po, o.s = 0.1, and 
i? = 10 km, we estimate the viscous timescale for a SS as 



77 



5.3 X 10^ Tg^/^ s 



(3) 



This is characteristic timescale of the angular momentum 
transport within a strange star. 



Local compression or decompression of the SQM drives 
it from the flavor equilibrium, because chemical poten- 
tials for the ultrarelativistic up and down quarks have 
weaker density dependence than the chemical potential of 
the massive strange quarks. Therefore, a baryon density 
perturbation implies /is 7^ /i^. The non-equilibrium reac- 
tions, changing the concentrations of the down and strange 
quarks, and driving the system toward equilibrium, are 
irreversible, and therefore imply dissipation. The dissipa- 
tion rate depends sensitively on the strange quark mass. 
Of course, dissipation vanishes in the (unphysical) limit of 
TOs — > 0, where the flavor symmetry is enforced. The dis- 
sipation rate in a pulsating SQM depends on the tempera- 
ture, density, frequency of pulsations w, and on the density 
oscillation amplitude 5n\j. The time- averaged dissipation 
rate can be represented by the (effective) bulk viscosity of 
SQM. The formula for C was derived by Madsen (1992), 
and can be written as 



C = 9 X 10- 
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V 100 MeV J VlO^s 

1\21 



\ 4 



P_ 

PO 



X [l + bTg^Uj/lO^ S-1) 



(5) 



where the dimensionless coefficient b is proportional to 
the coefficient (3 of Madsen (1992), and the terms depend- 
ing on Sn^ were assumed to be unimportant (see also 
Andersson, Jones, Kokkotas, 2002). The second term in 
square brackets is important for T <; 10^" K, where it 
leads to the quenching of opening in this way the possi- 
bility of the r-mode instability in very hot rotating strange 
stars (Madsen 2000). 

The terms depending on Sn^ become important when 



> 2 



(6) 



Such a condition can be realized at low temperatures 
T < 10^ K only provided 5n\,/n\, in the r-mode is siz- 
able, 5n\^/n\^ 0.1: this could happen only for a large- 
amplitude r-mode instability. In what follows, the second 
term in the square brackets will be omitted. 

Consider now T < 5 x 10® K, when the square bracket 
in Eq. (|SJ| can be replaced by one. Using results of Madsen 
(2000) and assuming a mean density of strange star ipo, 
we get the following rough estimate for the viscous damp- 
ing timescale of the r-mode, 



T < 5 X 10 -" K 



nisC 



100 MeV 



/ 



1000 Hz 



T, 



9 s. (7) 



The damping timescale depends strongly on , which is 
a renormalized strange quark mass (see Sect. 3.1). 



6.2. Viscous timescales connected with bulk viscosity 6.3. Cooling timescale 

Bulk viscosity C of the SQM is due to the non-equilibrium The SS cools rapidly, with 

T cvT 



non-leptonic strangeness changing process 
u + d ^ s + u . 



(4) 



'^cool 



(8) 
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where cy is the specific heat of SQM (at a fixed volume), 
and Qi, is the neutrino emissivity. Eq. © can be easily 
integrated to give the time (age) needed to cool SS to 
Tg = T/10^ K, 

_ 108 

icool — S . (9) 



6.4. Damping of the Jacobi-like instability 

Lindblom & Detweiler (1977) have shown that if both vis- 
cosity and gravitational radiation reactions are taken into 
account, rotating stars can be stable beyond the Jacobi 
bifurcation point located in Sec. This stability can be 
seen as resulting from a competition between the CPS 
Dedekind-like mode and the Jacobi-like mode. It is gov- 
erned by the ratio (deduced from Eq. (8) of Lindblom & 
Detweiler 1977) 



X = 4.41x10^^^ (- 



R 



10 km 



M 



1.4 M, 



,(10) 



where v ~ rj/ p is the kinematic shear viscosity, li X ^ I 
the star is unstable with respect to the Jacobi-like mode 
(viscosity-driven instability), whereas if X <C 1, it is un- 
stable with respect to the Dedckind-like mode (CPS in- 
stability). In the intermediate regime, X ^ 1, the star 
is stable well beyond the bifurcation point. The formula 
(|10|) has been established for incompressible stars. Taking 
into account the compressibility does not change it sig- 
nificantly, since for polytropes of adiabatic index 7 this 
merely adds a factor (7 — 6/5)^/(7 — 1)^ in front of it (Lai 
& Shapiro 1995). 

Putting in Eq. (|10|) the value of ly deduced from the 
shear viscosity given by Eq. ||2J) leads to 



X = 2 X 10" 



R 



0.1 

a 
2 



10 km 



5/9 



10^ K 
T 



5/3 



1.4 M, 



(11) 



The fact that X 1 for astrophysicaly interesting tem- 
peratures T> lO^K (see eq. El means that the viscosity- 
driven instability cannot develop in rotating strange stars 
described by the MIT bag model, even if a star rotates 
rapidly with T/\W\ > iT/\W\)crit. 

6.5. The instability windows for r- modes 

The CPS-instability of r-modes in strange stars has been 
studied by Madsen (2000) and Andersson et al. (2002) 
in the Newtonian limit. They found the critical rotation 
frequency for a given stellar model as a function of tem- 
perature solving the following equation 

- + - + - = 0, (12) 

Tr 



Where r. ^ _22(47) (d^) " (^)-' (^X ^ 
(Kokkotas & Stergioulas 1999) is the r-mode growth time- 



scale, taking into account only the leading current mul- 
tipole and assuming a constant density profile (or n=l 
polytrop). According to this study, the instability oper- 
ates in a very different temperature window than for NS 
(see Pig 1 in Madsen 2000 for the dependence of critical 
spin frequency on internal temperature for the MIT bag 
model strange star with tUsC^ = 100 and 200 MeV). Since 
for SS bulk viscosity becomes weak at both very low and 
very high temperatures, the r-mode instability could oc- 
cur in temperature windows T > 10^" K (4 10^ K) and 
T < 10^ K (2 10^ K) for nisC^ = 200 MeV, shown as the 
hatched area in our Pig. 9 (for mgC^ = 100 MeV above 
dotted line) , if the cooling time is long enough. Andersson 
et al. (2002) suggested that the r-modes in a strange star 
never grow to large amplitudes and emit a persistent grav- 
itational wave signal. 

6.6. Astrophysical scenarios 

Let us consider two evolutionary scenarios for rotating 
strange stars. The first scenario refers to a newly born, 
rapidly rotating SS, cooling due to neutrino emission. The 
second scenario is that of a SS in a close binary system, 
which heats and spins-up due to accretion of matter. In 
both cases we assume that strange quark matter (SQM) 
is normal (non-superconducting). 

There are important differences between the two as- 
trophysical scenarios mentioned above. The first one is 
that the energy radiated by a newborn relativistic star is 
supplied by its rotational kinetic energy, if it acquired a 
sufficiently large amount of rotational energy during its 
formation. Por an accreting compact star it is supplied 
by the accreted matter. An old strange star in a close bi- 
nary system, accreting matter from its companion, may 
be spun up until the ratio of is high enough to 

allow for the symmetry breaking. A steady state regime 
is achieved when the total accreted angular momentum 
is evacuated via gravitational radiation (Wagoner 1984, 
Papaloizou & Pringle 1978), this process being known as 
forced gravitational emission or the Wagoner mechanism. 
In this case the frequency of the emitted gravitational 
waves depends upon the type of triaxial instability. Por 
the main I — m = 2 r-mode the frequency of the emerging 
gravitational waves is (see e.g. Stergioulas 1998) 

/•r— mode ^ 



J g 



Por the Jacobi-like bar mode instability, the frequency 
of gravitational waves is simply twice the rotational fre- 
quency: 

j?visc o J" 

J gw .1 ' 

If, in addition to the measure of /g^ by a gravitational 
wave detector, one knows / from electromagnetic obser- 
vations (e.g. X-ray pulsations), one can determine the kind 
of secular rotational instabilities acting on the star, and 
consequently have some insight into the value of the physi- 
cal parameters (viscosity, temperature) inside the compact 
object. 
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6.6.1. Newly born rapidly rotating strange star 

A SS is born from a proto neutron star with an initial 
temperature of ^ 10^^ K, and cools via neutrino emission 
(shown in a schematic way as the N track in Fig. 9). Of 
two basic types of SQM viscosity, the bulk one is crucial 
for stabilizing rapid rotation with respect to the r-mode 
instability (which develops on a tiniescale r^.). The region 
of the r-mode instability in the / — T plane strongly de- 
pends on the value of mg, because cx m~^. Had we taken 
TOgC^ = 100 MeV, the r-mode instability of newly born SS 
would extend to much lower T and /. We use renormalized 
s quark mass nisC^ = 200 MeV consistent with particle 
data tables (see Sect. 3.1). Then the r-mode is unstable 
within the right hatched strip in Fig. IHl Due to the strong 
temperature dependence of r<^, stability of rotation is re- 
gained as soon as T decreases below 10^" K, which takes 
only a few seconds (see eq.lH)). A rapidly cooling young hot 
strange star would not be spun down significantly passing 
the right r-mode instability window. 

At T = IQS K, t,, ~ 10^ s < tcooi(10^ K) = IQS s 
and a star could rotate rigidly. Since A" << 1 (due to a 
very short gravitational damping time-scale) the viscosity 
driven instability cannot develop even if the newly born SS 
left the r-mode instability strip with r/|VK| > {T/\W\)cnt- 

Several years after its birth a star will reach the low- 
T (left) (when temperature is a few times 10^ — 10^ K) 
r-mode instability strip. During its cooling it will lose the 
excess of angular momentum via GW radiation until it 
leaves the low-T r-mode instability window. 



6.6.2. Spin evolution of an old accreting strange star 

Consider an old SS which starts to accrete matter in a 
close binary system. Its initial internal temperature was 
very low, but energy release associated with accretion will 
induce heating of SS interior. To be specific, we consider 
accretion of a hydrogen-rich plasma. There are two main 
channels of heating by accretion. The first one is connected 
to plasma falling onto a bar SS surface, and the second one 
involves formation of the crust and thermonuclear burning 
of compressed plasma (Haensel & Zdunik 1991). In both 
cases, additional heat release occurs at the SQM core edge, 
where absorption and deconfincment of protons is associ- 
ated with exothermic irreversible reactions p — > 2u + d, 
u + d — > s + u. The energy release per one deconfined 
proton is Qs- In what follows we will assume qs — 50 MeV. 
After some 10^ yrs SS, accreting matter at a constant 
rate M, it will reach the steady thermal state, in which 
heating due to accretion and deconfinement is balanced 
with neutrino losses from the SS core (see Fujimoto et al. 
1984). Notice that at a fixed accretion rate (which will 
be expressed in the units of lO^'* Mq yr~^), Af_8, reach- 
ing of the steady state corresponds to accretion of only 
10~3 M_8 Mq (H-track in Fig. Equating heat gams 



and losses, one can estimate internal temperature of ac- 
creting SS (Haensel & Zdunik 1991), 

T = 1.1 X 10® (^CM^M-sY^^ Rq^^"^ K , (13) 

where M^. = M/Mq and Rq = i?/10^ cm. The parame- 
ter ^ varies from ^ — 0.2 in the limiting case when the 
heating results exclusively from the deconfinement pro- 
cess, to ^ = 1.2 when all gravitational energy release by 
infalling protons flows into the SS core (Haensel & Zdunik 
1991). Putting M, = 1.6 and ijg = 1, ^ = 0.2, we get 
T = 0.9 X 10® K. 

Further accretion of some (0.3 — 0.4) Mq during some 
5 X 10^ yrs at AI-s — 1 could spin-up the SS to frequency 
/ ^ IkHz, at a nearly constant temperature (Zdunik et al. 
2002) (the vertical S-track in Fig. EJ. Along the S-track, 
for TTis — 200MeV, a star is stable with respect to the r- 
mode instability because of its huge bulk viscosity. It can 
be accelerated to high rotational frequency reaching the 
critical frequency for viscosity driven instability (short- 
long dashed fine in Fig. EJ. After crossing / = /dit = 
(0.9 — 1.2) kHz (the precise value depending on the SS 
mass and the detailed form of the SQM EOS) the SS is 
stable to the bar mode formation since the characteristic 
growth time-scale of viscosity driven instability ~ a 
month >> than tens of seconds - the gravitational wave 
damping time. In the case of rus = 200MeV a strange star 
can be accelerated to very high rotational frequencies until 
it reaches one of the upper limits (low-T r-mode instability 
region, dynamical instability limit, mass-shedding limit). 

Let us stress the importance of the value of mg for 
such an astrophysical scenario. For high mass of strange 
quark > 200 MeV a star can be accelerated to very high 
frequencies without emitting gravitational waves. If on the 
contrary the strange quark mass is mg < 100 MeV the low- 
T r-mode instability window (above the dotted line in Fig. 
inj, it can be reached at a rotation rate of 0.4-0.6 kHz. The 
rotating SS is expected to be in a quasi-stationary state, 
in which the angular momentum losses via gravitational 
radiation are balanced by the accretion torque. In this 
steady state an accreting SS is an emitter of a persistent 
gravitational wave signal. 

6.6.3. The role of the magnetic field 

Our results could be changed if we take the mag- 
netic field into account. The presence of a strong mag- 
netic field in SQM could also contribute to the an- 
gular momentum transport and therefore influence dif- 
ferent kinds of gravitational instabilities. The electrical 
conductivity of the SQM is huge. Using the formula 
derived by Heiselberg & Pethick (1993) we get a ~ 
(0.1/a)5/3 {p/pof^ Tg^'^ 10^1 s-i. Therefore, the mag- 
netic field is frozen in SQM which combined with the pres- 
ence of rapid and differential rotation could lead to addi- 
tional viscous mechanisms due to the magnetorotational 
instability. 
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Fig. 9. Regions of the r-modc instability (hatched region 
with r-label) for a canonical mass strange star assuming 
msC? = 200 MeV in the rotational frequency- temperature 
plane. The r-mode instability domain for rrigC^ = 100 MeV 
is shown for comparison (regions above dotted lines). The 
r-mode instability windows are taken from Madsen 2000. 
The short-long dashed horizontal line corresponds to the 
critical rotation rate for viscosity driven instability assum- 
ing canonical stellar mass, mgC^ = 200 MeV and low MIT 
bag constant 40MeV/fm'^. The H-S tract corresponds to 
the evolution of an old accreting SS and N tract to a newly 
born strange star. 



Although a high conductivity is a key element of com- 
pact stars (neutron stars, strange stars) physics there have 
been only a few discussions of the effect of magnetic fields 
on triaxial instabilities (Bonazzola et al. 1996, Spruit 1999, 
Rezzolla et al. 2000, Ho & Lai 2000). It was suggested by 
Bonazzola, Frieben & Gourgoulhon (1996) that the mag- 
netic field can prevent the Dedekind-like mode. Assuming 
infinite conductivity, shear of the fluid from differential 
angular velocity generates some toroidal magnetic field in 
addition to the poloidal one (Bonazzola & Marck 1994). 
The excess of kinetic energy contained in differential rota- 
tion with respect to rigid rotation can be efficiently con- 
verted into magnetic energy, enforcing rigid rotation. 

Ho & Lai (2000) have studied the combined effects of 
magnetic braking and gravitational radiation generated by 
r-mode instability on the spin evolution of young rapidly 
rotating neutron stars. They suggested that a strong mag- 
netic field can drive the instability just like gravitational 
radiation reaction. The question of whether an unstable r- 
mode leads to differential rotation, and whether the mode 
can be prevented from growing by the magnetic field, was 
discussed by Rezzolla, Lamb and Shapiro (2000). They 
studied how so-called Stokers drift affects the magnetic 
field of the star, and what the backreaction on the mode 
may be. It was estimated that the instability could op- 



erate in young neutron stars and recycled ones provided 
that they spin fast enough. 

7. Summary 

We have performed numerical investigations of the Jacobi- 
like bar mode secular instability of rotating strange quark 
stars described by the MIT bag model. In particular, we 
have computed the domain of the parameter space of sta- 
tionary rigidly rotating configurations in which triaxial 
instability develops. It has been found that, contrary to 
neutron stars, the triaxial instability can set in quite far 
towards the mass-shed limit and is allowed for small mass 
stars as well if the viscosity is high enough to damp out 
any deviation from rigid rotation. The lower mass is the 
lower rotational frequency at which Jacobi-like instability 
can develop. The viscosity driven instability does not set 
an upper limit on the spin rate (but on r/|VK| or oblate- 
ness) of strange stars - the configuration marginally stable 
against non-axisymmetric perturbation can rotate with a 
frequency higher than mass-shedding configurations (the 
maximal rotating configuration is not a Keplerian one 
for SS). A normal and low mass supramassive strange 
star gaining angular momentum always slows down before 
reaching the Keplerian limit. For a high-mass supramas- 
sive strange star sequence, the mass-shedding configura- 
tion is the one with the lowest rotational frequency in the 
sequence. The instability, if relevant, can impose a strong 
constrain on the lower limit of the frequency at the inner- 
most stable circular orbit around rapidly rotating strange 
stars. The results arc robust for all linear self-bound equa- 
tions of state. The critical rotation frequency for viscosity 
driven instability in the case of strange stars described by 
the MIT bag model is high, 0.8 — 1.2 kHz (depending on 
the stellar mass and the detailed form of the SQM EOS). 

We have discussed astrophysical scenarios when triax- 
ial instabilities (r-mode and viscosity driven instability) 
could be relevant in strange stars described by the stan- 
dard MIT bag model of normal quark matter (any kind of 
superconductivity being excluded). We have shown, tak- 
ing into account actual values of shear viscosity, that vis- 
cosity driven instability cannot develop in any astrophys- 
icaly relevant temperature windows. 

Whether a strange star can be a good emitter of grav- 
itational waves or not depends on CFS instability. In our 
discussion we focused on the most promising CFS insta- 
bility (r-mode). The spin evolution of strange stars and 
emission of gravitational waves strongly depend on the 
mass of the strange quark (caused by the strong depen- 
dence of bulk viscosity on the mass of a strange quark). If 
we take the strange quark mass to be rrisC? ~ 200 MeV or 
higher (consistent with renormalization and particle data 
tables) then r-mode instability is relevant only in a nar- 
row temperature window and seems to be unimportant in 
Low Mass X-Ray Binaries. In this case an old star can be 
spun up in LMXBs to very high rotation rates and a low 
frequency at the marginally stable orbit can be obtained 
(~ IkHz). 
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If on the contrary the strange quark mass is low, 
TTigC^ ~ 100 MeV, then a star accelerated in LMXBs to a 
frequency of 0.4 — 0.6 kHz can be a good emitter of persis- 
tent gravitational waves (see Andersson, Jones, Kokkotas 
2002 for a discussion) that ought to be detectable with 
large-scale interferometers. In this case the r-mode insta- 
bility can impose strong constraints on the lower limits 
on rotational period and then on the maximal orbital fre- 
quency in LMXBs. The lowest frequency at the marginally 
stable orbit could be obtained for slowly-rotating massive 
strange stars. 

A young strange star would not be spun down signifi- 
cantly by the r-modes during its first few months of exis- 
tence (because of its very rapid cooling rate) until the star 
has cooled to a temperature of a few times 10^ — 10® K. 

The astrophysical scenario will change if values of vis- 
cosity change or/and different physical mechanisms e.g. 
magnetic field will be taken into account. 
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